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Abstract 

We consider a two flavor Polyakov-Nambu- Jona-Lasinio (PNJL) model where the Lagrangian 
includes an interaction term that explicitly breaks the Ua(1) anomaly. At finite temperature, 
the restoration of chiral and axial symmetries, signaled by the behavior of several observables, is 
investigated. We compare the effects of two regularizations at finite temperature, one of them, 
that allows high momentum quarks states, leading to the full recovery of chiral symmetry. From 
the analysis of the behavior of the topological susceptibility and of the mesonic masses of the axial 
partners, it is found in the SU(2) model that, unlike the SU(3) results, the recovery of the axial 
symmetry is not a consequence of the full recovery of the chiral symmetry. Thus, one needs to 
use an additional idea, by means of a temperature dependence of the anomaly coefficient, that 
simulates instanton suppression effects. 
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I. INTRODUCTION 



Spontaneous breaking and restoration of chiral symmetry is one of the most fundamental 
aspects inferred from quantum chromodynamic theory (QCD) and is manifested through 
the equation of state of hot and dense hadronic matter. The Nambu-Jona-Lasinio (NJL) 
model is one of the most studied effective models of QCD at finite temperature and has been 
applied successfully to reproduce the low-energy phenomena, although this model does not 
support some features of QCD, such as quark confinement. 

With the intention of broadening the applicability of effective QCD models, in particular, 
allowing for the implementation of additional features of the QCD phenomenology, the 
Polyakov Loop is considered within the NJL model, that so is endowed with a mechanism 
that allows the analysis of the expected confinement /deconfinement phase transition in the 
hadron-quark system. This extended model, the so called Polyakov-Nambu-Jona-Lasinio 
(PNJL) model, first implemented in Ref. l|, provides a simple framework which considers 
the chiral and the confinement order parameters. In fact, the NJL model describes the 
interactions between constituent quarks, hence providing the correct chiral properties. The 
static gluonic degrees of freedom, introduced in the NJL Lagrangian through an effective 
gluon potential in terms of the Polyakov loop, take into account features of deconfinement 
lH4j. A good approach to reproduce lattice results is given by the coupling of the quarks 
to the Polyakov loop. This leads to a reduction weight of the quark degrees of freedom 
as the critical temperature is approached from above (interpreted as a manifestation of 
confinement). 

Besides the chiral symmetry, another important symmetry, explicitly broken in the QCD 
Lagrangian, is the axial Ua(1) symmetry Q], which might also be restored at extreme temper- 
atures. In fact, there are indications from the lattice calculations that, at high temperatures, 
effects arising from the Ua(1) breaking are strongly suppressed. This suggests an effective 
restoration of the Ua(1) symmetry. An important issue that the present paper intends to 
analyze is the possible restoration of axial and chiral symmetries, which observables may 
give indications regarding this restoration, and whether there is an interplay between both 
restorations or not. By this we will use the PNJL model in the SU(2) sector but with the 
inclusion of a term on the Lagrangian that explicitly breaks the Ua(1) symmetry. It is worth 
noting that the effects of the anomaly have been intensively studied in the NJL model in 



SU(3) with the 't Hooft term 6Ml0l|. and more recently, the possibility of axial symmetry 
restoration and its effects on diverse observables jf], 9K11]. The model in SU(2), due to 
its simplicity, allows to isolate some aspects of the problem, reason why its study can lead 
to a relevant contribution for the understanding of physics associated to the breaking and 
restoration of the Ua(1) symmetry. Interest in studying this issue in two flavor models is 
manifested in recent works [l^ . 

The same model has been used in Ref. 13J which particularly looks into the nature of the 



superfluid phase. Here we will focus on the topological susceptibility, the chiral and axial 
partners and the comparison between the results of SU(2) and SU(3) models. 

We will investigate the behavior of the following observables with temperature: quark 
condensates, topological susceptibility and it, a, a® and r\ mesons. Two approaches will be 
followed in order to understand the mechanism of restoration of axial symmetry: we will 
consider different degrees of Ua(1) symmetry breaking in the vacuum, by using different 
values for the strength of the anomaly that are kept constant with varying temperature; 
alternatively, we fix the degree of symmetry breaking in the vacuum and allow the coupling 
strength to decrease with temperature. We will also consider the effects upon the observables 
behavior of using two types of regularization: 

- Regularization I: as the cutoff A is only necessary to regularize some integrals in the 
vacuum and it is not necessary at finite T, it is considered infinite as soon as thermal effects 
are considered (A — > 00). This has the effect of allowing for high quark moments; 

- Regularization II: the cutoff is taken always with a finite value, at zero or finite T 
(A = const.), and all the momenta couple with the same strength up to a cutoff momentum 
A ~ 1 GeV. 

In Sec. II we will introduce the formalism for a two flavor NJL type model with two 
different interacting parts that allow to disentangle chiral and axial symmetries. The gap 
equations and the mesonic propagators will also be discussed, as well as the breaking of 
the axial symmetry and the topological susceptibility. The Polyakov loop extension of the 
model will be explained as well. The different types of regularization will be presented in 
the context of the model extension for finite temperature. We redirect the details of the 



calculations to the appendices included in a previous paper [10|. In Sec. Ill we present 
and discuss our results for the phase transition and also for the temperature behavior of 
the meson masses and the topological susceptibility, aiming at discussing restoration of 



symmetries. We conclude in Sec. IV with a brief summary. 



II. FORMALISM 

The aim of this section is to present the mathematical formalism of the N JL model and its 
extension to the PNJL model. In the present work, a two flavor, three color, quark system 
will be studied, and it will be explicitly introduced a term on the Lagrangian that acts as 
a source of anomaly. The anomaly of the SU(2) model is already present in the original 
NJL model Q, flol through a Fierz transformation and has been considered in other works 
as well (see Ref. [16| and Refs. there in), the strength of the anomaly coupling being the 
same as the chiral coupling. In the present work we have an extra term that ensures an 
independent mechanism of axial symmetry breaking. 



A. A SU(2) Nambu Jona-Lasino type model 

The Lagrangian that will be used here is a SU(2) version of the NJL model taking into 
account an additional term that, although being a chiral invariant, explicitly breaks the 
axial symmetry. 

We start with a NJL type model with two flavors defined by the following Lagrangian: 



£ = q(i$ - m)q + Ci + C 2 , 



(1) 



with two different interacting parts 



£1 = 9i (QQ) + ioilsTq) + (qfq) + (qi^qY 



(2) 



£2 = 92 (qq) + (qil5Tq) - {qfq) - (qi^qY 



(3) 



The quark fields q = (u, d) are defined in Dirac and color fields, respectively with two 
flavors, Nf = 2 and three colors, N c = 3, the coupling constants g\ and g 2 have dimension 
energy -2 , and rh = dia,g(m u ,m ( i) is the current quark mass matrix. 

The Lagrangian is chiral invariant in the limit where the current quark masses vanish. 
Both terms Ci and C 2 are invariant upon SU(2)i®SU(2) fl (8)U(l) type transformations, but 



the £2 component makes the Lagrangian non-covariant upon U^(l) transformations. The 
L 2 term, that may be represented in the form of a determinant: 



£2 = 2g 2 det [q(l + y 5 )q] + det [q(l - 75)?] 



(4) 



can be identified as an interaction induced by instantons, according to 't Hooft, and it 
explicitly breaks the axial symmetry even in the chiral limit. 

In this model, while the vector current is conserved, the axial isovector current is only 
conserved in the chiral limit and, as a consequence of the £2 term in the Lagrangian, the 
isoscalar axial current is not conserved even in the chiral limit. In fact: 



= 2m(g 2754) + 8g 2 {qrq){qi^fq) - (qq)(qi^ 5 q) 



The last equation is the equivalent of the QCD 4-divergence of j§\ 



d^ = 2m(q i l5 q) + 2N F Q(x), 
which allows to identify the topological charge Q(x) in our model as: 



(5) 



(6) 



Q(x) = 2g 2 (qrq)(qi-f 5 fq) - (qq)(qijsq) 



(7) 



or, equivalently, as 



Q(x) = 2g 2 det [(q(l - 7s )g - det [(q(l + 75) 



(8) 



This equation shows the important role played by the coupling constant g 2 in our analysis 
bearing in mind that the physical effects of the Ua(1) anomaly are only manifested with 
non-zero topological charge. 

Having identified the topological charge, the topological susceptibility may be calculated 



llj], as it will be shown in Sec. II. C. 



The four fields in the original Lagrangian can be rearranged in the form: 



£i + £ 2 = | 



(qq) + {qi^fqY 



9a 

2 



{qfqf + (qi^qf 



(9) 



where g s = 2 (5-1 + g 2 ) and g a = 2 (g x - g 2 ). 

From the last expression it is easy to understand the reason why the Lagrangian supports 
four mesonic channels: O a — 1, O n — ij^T, O r] = 275 and = lr. After integrating the 



generating functional over quark fields and using a, n, rj, do as auxiliary fields, the following 
effective action is obtained: 



l eff 



( -, -> \ a + 71" V + Oo 

-iTi In - id - m + a + i^m + 27577 + ra . (10) 

V / 2q s 2q n 



2g s 2g a 



A standard calculation leads straightforwardly to the gap equation and to the meson 
propagators. The interaction term of the Lagrangian with g s coupling is responsible for the 
calculation of the propagators of the chiral partners (ir , a), while the term associated to g a 
allows for the chiral partners (77,00). Since we consider equal quark masses, flavor mixing 
effects induced by the axial symmetry breaking are not visible. 

The following gap equations are obtained: 

M l =m i -2g s (qq) i , (11) 

where one identifies % = u,d and Mj as the constituent quark mass. The quark condensates 
are determined by 

(qq)i = ~^Tr 1 = -zTr^(p), (12) 

being Si(p) = (j/ — Mj + ie)~ l the propagator of quarks. 

The mass spectra of the mesons is obtained by the analysis of the pole structure of the 
meson propagator, given by 

l-2g s , a Tl M (q 2 = M 2 M )=0. (13) 

where 

n M (q2) = i f ^Tr [O m S(p + q) O m S(p)} (14) 

is the polarization operator for the quark-antiquark system regarding the channel with quan- 
tum numbers {M.} in the mesonic sector. As mentioned above, g s is related to tt and a 
mesons and g a to rj and Oq mesons. 



B. Extension to the PNJL model 

Now we include the Polyakov loop and its effective potential to the NJL type model 
described above. The Lagrangian of this SU(2)®SU(2) quark model with explicit chiral 
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a ai a 2 63 
3.51 -2.47 15.2 -1.75 

TABLE I. Parameters for the effective potential in the pure gauge sector. 



symmetry breaking where the quarks couple to a (spatially constant) temporal background 
gauge field (represented in term of Polyakov loops) is given by 0, [13] : 

Cpnjl = q(i^ D a - m)q+^ \(qq) 2 + (qi^rq) 2 } + y \{qfq) 2 + [qi^qf] 

-U($[A],$[A];T). (15) 

The quarks are coupled to the gauge sector via the covariant derivative D 11 = <9 M — iA^. 
The strong coupling constant g strong has been absorbed in the definition of A^: A^(x) = 
9 strong A^(x)^f- where A% is the SU C (3) gauge field and A a are the Gell-Mann matrices. 
Besides in the Polyakov gauge and at finite temperature A^ = 5qA° = —id^A 4 . 
The Polyakov loop $ (the order parameter of Z3 symmetric/broken phase transition in pure 
gauge) is the trace of the Polyakov line defined by: $ = -^((Vexpi dr A 4 (x,t) ))„. 

The pure gauge sector is described by an effective potential U ($[A], $[A\; T) chosen to 

reproduce at the mean-field level the results obtained in lattice calculations: 

U f$ T) a(T) - 

v y4 1 ; = 1^$$ + 6(T)ln[l - 6$$ + 4($ 3 + $ 3 ) - 3($$) 2 ], (16) 

where 

a (T) = a + ai fif) + a 2 (^j and b(T) = h (^j ■ (17) 

The effective potential exhibits the feature of a phase transition from color confinement 
(T < Tq, the minimum of the effective potential being at $ = 0) to color deconfinement 
(T > T , the minima of the effective potential occurring at $ 7^ 0). 

The parameters of the effective potential U are given in Table [B These parameters have 
been fixed in order to reproduce the lattice data for the expectation value of the Polyakov 



19] 



loop and QCD thermodynamics in the pure gauge sector 

The parameter T is the critical temperature for the deconfinement phase transition 
within a pure gauge approach: it was fixed to 270 MeV, according to lattice findings. This 
choice ensures an almost exact coincidence between chiral crossover and deconfinement at 
zero chemical potential, as observed in lattice calculations. 

The PNJL grand canonical potential density in the SU/(2) sector can be written as {^H]: 



$, M; T, //) = W ($, $, T) + 2g s N f (q^) 2 - 2N C N f J E p 

- 2N f T J ^3 (z+iEk) + , (18) 

where E { is the quasi-particle energy for the quark v. Ei — \/p 2 + Mf, and z J and z§ are 
the partition function densities. 

The explicit expression of z% and are given by: 

z+(Ei) = Tr c ln [l + lV^+m)] = hi {l + 3 (<5 + $ e -/W*0) e -/^+*0 + e * + ")}, 

(19) 

2 $ = Tr c ln [1 + Le-K Ei -ri] = In { 1 + 3 ($ + $ e -^^-^) e^*^ + e - 3 ^-M) } . 

(20) 

A word is in order to describe the role of the Polyakov loop in the present model. Almost 
all physical consequences of the coupling of quarks to the background gauge field stem from 
the fact that in the expression of z$, $ or $ appear only as a factor of the one- or two- 
quarks (or antiquarks) Boltzmann factor, for example e~^ Ei ~^ and e~ 2 ^ Ei ~v) . Hence when 
$, $ — > (signaling what we designate as the "confined phase") only e -3 ^^* - ' 4 ' remains 
in the expression of the grand canonical potential, leading to a thermal bath with a small 
quark density. At the contrary $, $ — > 1 (in the "deconfined phase") gives a thermal bath 
with all 1-, 2- and 3-particle contributions and a significant quark density. 

This formalism, presented here for completeness in the grand canonical approach, will be 
employed in the present work with \i = 0. This condition implies $ = $ = 0. 



C. The topological susceptibility 

The topological susceptibility, x, is an essential parameter for the study of the breaking 
and restoration of the Ua(1) symmetry. 

The topological susceptibility is defined as: 

X = j d 4 x (0\TQ(x)Q(0)\0) c , (21) 

where c means connected diagrams and T the time order operator. Using the definition of 
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Eq. ([7]) it can be seen that the operator Q(x)Q(0) in this model reads: 



Q(x)Q(0) = Agl [q(x)q(x)][q(x)j 5 q(x)] - [q(x)rq(x)][q(x)^Tq(x)} 
[g(0)g(0)][g(0)7 5 ?(0)] - [g(0)fg(0)][g(0) 75 rg(0)] 
= 4(G«) 2 [{q(x)q(x)} [q(x) l5 q(x)] [g(0)g(0)] [3(0)759(0)] 

- [q(x)q(x)] [q(x) l5 q(x)} [g(0)f g(0)] [g(0) T5 f g(0)] 

- [q{x)fq{x)} [q(x) l5 rq(x)} [g(0)g(0)] [g(0) T5 g(0)] 

+ [q{x)rq{x)} [q{x) lh Tq{x)\ [g(0)fg(0)] [g(0) 75 fg(0)] 



(22) 



Taking into account only the connected diagrams of order 1/N C (see Fig. [T]) and fol 
lowing a similar approach to 
susceptibility: 



11] . we arrive at the following expression for the topological 

4/i 



X = 4N f g 2 2 (qq) 



(23) 



1 - SgJi ' 

As already mentioned, one of the aims of this paper is to perform a comparison between 
the behavior with temperature of the topological susceptibility in the framework of the 
SU(2) and SU(3) PNJL models with explicit axial symmetry breaking. As we will see, some 
meaningful differences will appear, so it is useful, in order to understand these differences, 
to compare both expressions. For the sake of simplicity we will write the SU(3) expression 
for the case of equal quark masses m u = = m s : 



X = ±N f (g D (qq)Y (qq) 2 - 



4/i 



8 -Poo h 

where -Pqo^. ids — %9d (qq)) (f° r details concerning the SU(3) model see the Appendix 



(24) 



and Ref. 



2 ll|). A comparison shows that the expressions for quantities with effects of the 



anomaly in SU(2) depend on a constant, 5% while the corresponding expressions in SU(3) 
depend on a constant multiplied by a quark condensate, go{qq)- The same replacement 
holds for the gap equations and for the propagators of the chiral partners (r], ao). 

We will return to this point again when discussing the behavior of several observables with 
temperature where it will be shown that some effects are different in SU(2) and SU(3) due 
to the fact that in the last model the "effective" constant gD e ff = 9o(qq) already depends 
on temperature through the quark condensate, while its corresponding in SU(2), g 2 , remains 
constant, unless we explicitly allow a variation with temperature. 
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c) 




FIG. 1. Connected diagrams for the loop ((a) and (b)) and for the ring contributions (c). Only 
(b) and (c) are of order 1/N C . 

III. RESULTS 

A. Vacuum properties and parameters fixing 

The present PNJL model has four parameters in the NJL sector: m, A, gi and gr 2 - We 
choose to adjust the parameters in vacuum by fitting to well known experimental data 
or lattice values: the mass of the pion, its decay constant, the quark condensate and the 
topological susceptibility, x-, that are shown in Table [III The masses of the a, r\ and a 
mesons come as outputs. 

Our main purpose is to study the behavior of several observables with temperature, 
in particular those that might signal the chiral and deconfinement phase transitions, and 
the effective restoration of chiral and axial symmetries. For this we will use two different 
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u 

[MeV] 


(gg)V3 
[MeV] 


m w 
[MeV] 


m a 
[MeV] 


m v 
[MeV] 


m ao 
[MeV] 


x 1/4 

[MeV] 


Model 


93 


-300 


140.2 


803.7 


704.5 


919.8 


180.8 


Exp. /Latt. 


92.4 


-270 


135.0 


400-1200 


547.3 


984.7 


180 



TABLE II. Numerical values for the calculated observables: f n , (qq), the meson masses and the 
topological susceptibility, obtained with A = 590 MeV, GA 2 = 2.435, a = 0.2 and m = 6 MeV. 



approaches, that will be presented now in the vacuum, and its consequences for the behavior 
with temperature will be explored in the next subsection. The two scenario to be considered 
are: 

Scenario A - We will keep g\ and g 2 as independent parameters. At finite temperature 
we may allow g 2 to be temperature dependent but g\ is kept constant. This is equivalent to 
the usual treatment of the SU(3) model; 

Scenario B - We redefine the coupling constants such as the set (^1,^2) will be replaced 
by (G, a) in the following parametrization: 



9i = G(l 



a) 



g 2 = Ga, 



(25) 



with a e {0, 1} 13|,|l6[. In this picture g\ and g 2 are not independent, but g s = 2(g 1 +g 2 ) = 2G 
will be kept always constant; on the contrary, g a = 2G(1 — 2a) varies with a. 

Consequently, when studying the effect of the degree of axial symmetry breaking, by 
choosing different values for a, only the topological susceptibility and the masses of 77, 
will be affected, since they depend explicitly on a; the n and a channels will not be affected. 

Let us now discuss Scenario B for the vacuum state considering special values for the 
parameter a: 

(1) a = 0: the coupling constant g 2 = and the axial symmetry is unbroken. 

(2) a — 1: the coupling constant g\ = and we have a pure instanton interaction, where 
Ua(1) is broken maximally. 



(3) a = 0.5: in this case g± = g 2 = G, the standard NJL model 16| is recovered, and only 
the a and n channels, are reproduced. 
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a 


X 1/4 (180) 


7?(547) 


a (984) 




[MeV] 


[MeV] 


[MeV] 








140.2 


803.7 


0.10 


150.7 


510.2 


856.7 


0.112 


155.5 


547.0 


863.7 


0.16 


170.6 


641.0 


893.2 


0.20 


180.8 


704.5 


919.8 


0.25 


191.5 


766.1 


956.3 


0.30 


200.7 


799.3 


996.4 


0.40 


216.0 


986.9 


1094.3 


0.50 


228.6 







TABLE III. Numerical values for the topological susceptibility, and masses of the mesons rj and a$ 
for different values of the anomaly parameter, a. 



In order to see the effects of the degree of Ua(1) symmetry breaking in the vacuum, we 
allow a to take values between and 0.5. The results are summarized in Table IIHI As 
it can be seen, for a = (no anomaly) the rj is degenerated in mass with the pion, the 
cio with a, and the topological susceptibility is zero, as expected. A small breaking of the 
axial symmetry (a = 0.1) is enough to break the degeneracy, raising the masses of rj (that 
comes close to almost its experimental value) and do, and to get a meaningful non- vanishing 
value of x 1 / 4 . Larger values of a yield larger values of x 1//4 and of the meson masses. When 
0.16 < a < 0.25 the value of x is within the range of values coming from lattice results 



22 



23]. 



1. Behavior with temperature 

Let us now discuss our results at finite temperature within the SU(2) PNJL model with 
anomaly. We aim at analyzing signals of possible restoration of chiral and axial symmetries 
and also to compare these results with those obtained within the corresponding SU(3) model 



21 



241 ] . First we analyze the characteristic temperatures for the chiral/deconfinement phase 



transition. An important question to address is whether 11^(1) is effectively restored at finite 
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FIG. 2. Comparison of the behavior with temperature of — (qq) 1 ^ 3 and the Polyakov field $ (left 
panel) and of the topological susceptibility (right panel) with 52 constant and two regularizations 
(/: A — > 00, //: A = const.) and with §2 a decreasing function of temperature and regularization 



II (a = 0.2). The lattice points for \ were taken from Ref. 



221. 



temperature and its possible relation with the restoration of chiral symmetry. To quantify 
this effect it is mandatory to analyze the behavior of the topological susceptibility, and of 
the mesonic excitations with temperature. Another important question in this concern is to 
discuss if x should be considered as an order parameter for the axial symmetry. A challenging 
question is whether the restoration of axial symmetry arises as consequence of the restoration 
of chiral symmetry, or if there is an independent mechanism for the suppression of instanton 
effects. To simulate the last effect, the anomaly coupling #2 may be assumed as a decreasing 
function of temperature. 

We start by discussing the behavior with temperature of observables that may signal 
restoration of chiral symmetry, deconfinement and restoration of axial symmetry in different 
situations. For this purpose, let us first analyze in Fig. |2] (left panel) the behavior of 
the quark condensate — (qq) 1 ^ 3 , of the Polyakov field $ (left panel) and of the topological 
susceptibility (right panel) with $2 constant and two regularizations (I: A — > 00, II: A = 
const.). From now on we will focus our discussion for the case a = 0.2 (the analysis is 
qualitatively similar for others values of a.) 
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or the Polyakov loop are similar to those 
, l^j], the critical temperature for the chiral 



The results for the quark condensate and_ 
within the SU(3) model. As it was shown in 
and deconfinement transitions, identified with the derivatives of the quark condensate and 
of the Polyakov field, respectively, become closer using regularization I. As shown in jjj, [2^ 
this regularization, together with the Polyakov loop, is very convenient to describe the 
thermodynamic properties of the system leading to a good agreement with lattice results. 
It has the disadvantage, for high temperatures, of leading to a sharp decrease of the quark 
condensate that, when using regularization /, vanishes and, unless an additional condition 
is implemented, becomes negative. 

In order to discuss the influence of an independent mechanism of suppression of instan- 
tons with temperature, we consider a third situation: the anomaly coupling is assumed to 
decrease with T, g2(T) = g 2 (Q)e~( T ~ 300 ^ 20 , and regularization II is used. Since Scenario A 
is considered, the variation of g 2 with T implies also a variation of g s = 2(g 1 + g 2 ), and there- 
fore the "mechanism of instanton suppression" influences the restoration of chiral symmetry 
and deconfinement. In fact, we can see from Fig. [2] that both — (qq) 1 ^ 3 and $ are affected 
by the decrease of g 2 . By examining the values of Table [IV] we see that, concerning the 
critical temperatures, both this case and regularization I with g 2 constant lead to a greater 
proximity of the deconfinement and chiral phase transitions as compared with regularization 
II and g 2 constant. We emphasize that the ansatz g 2 (T) leads even to equal values for T c x 
and T*. 

Observing the right panel of Fig. [2j we see that, as already shown in the SU(3) model 



2l|, PNJL calculation for the topological susceptibility nicely reproduces the first lattice 



points, a feature that is not verified in the NJL model [2lj. The effects of the introduction of 
high quark moments through regularization I is only relevant at high temperatures, allowing 
a decrease of x, that eventually becomes zero. The vanishing of x is & lso achieved with a 
decreasing with temperature of the strength of the anomaly (g 2 (T)). So, after analyzing 
both figures, it seems that two independent mechanisms, the allowance of hight momentum 
quark states and the instanton suppression, produce some similar effects: proximity between 
chiral and deconfinement phase transitions, and sharp decrease and vanishing of the topo- 
logical susceptibility. The question is whether the mentioned behavior of the topological 
susceptibility is enough to signal restoration of the axial symmetry. 

With rising temperature we expect the chiral SU(2)y® SXJ(2)a — SU(2)l®SU(2)^ sym- 
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regularization I 
(A -> oo) 


regularization II 
(A = const.) 


regularization II 
with g2(T) 


T c x [MeV] 


237 


287 


231 


Tf [MeV] 


219 


243 


231 


T c [MeV] 


228 


265 


231 


T* f [MeV] 


~ 260 


~ 350 


~ 245 


T e A ff [MeV] 






~ 285 



TABLE IV. Critical temperatures for chiral and deconfinement phase transitions (T* and Tf] T c is 
the average between both) and for the effective restoration of chiral and axial symmetries (T^j and 
T^jrj ), obtained for a = 0.2. The values in the two first rows were obtained with finite and infinite 
cutoff, respectively, and 32 constant; in the third row 52 is a decreasing function of temperature 
and the cutoff is kept constant. 

metry to be restored, and the chiral partners (tt, a) and (77, ao) become degenerate in mass, 
the effective restoration of chiral symmetry being then achieved. 

Concerning the axial symmetry, the vanishing of the mass splitting of the axial chiral 
partners (tt, ao) and (a, rj) is a necessary condition for the restoration of this symmetry. In 
order to discuss this problem, we are going to analyze the behavior with temperature of the 
topological susceptibility and meson masses by exploring two situations: (i) using the two 
regularizations already presented, but without allowing the anomaly coupling #2 to decrease 
with temperature, (ii) considering several degrees of violation of axial symmetry ab initio 
by using three different values for the anomalous parameter a. Once again we will consider 
a fixed a = 0.2 for T = 0. 

The temperature dependence of the mesonic observables is displayed in Fig. [3j Con- 
cerning the chiral partners (a, tt), we observe that in the vicinity of 260 (350) MeV for 
regularization I (regularization II) both masses become identical, a sign for the effective 
restoration of chiral symmetry. At a temperature slightly lower there occurs the conver- 
gence of the chiral partners (77, ao). The interesting result is that the restoration of the 
axial symmetry is not achieved: the masses of the two partners (it, a) and (rj, ao), although 
getting close at high temperatures, do not converge. Concerning the variation of a, the 
lower values of this parameter favor the restoration of axial symmetry, as expected, since 
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FIG. 3. Masses of the meson chiral partners (vr,cr) and (ij,ao) with varying temperature for three 
different values of the anomaly parameter a with two regularizations at finite T: regularization I 
(left panel) and regularization II (right panel). The dotted lines represent the qq threshold 2M U . 



the masses of the axial chiral mesons were already closer in the vacuum. Even when high 
momentum quarks are allowed (left panel) the masses of partners, although closer, do not 
converge (see Fig. [3j right panel), in spite of the vanishing of the topological susceptibility 
in this case (see Fig. HJ). 

This result is different from the one obtained in the SU(3) versions of NJL and PNJL 



models [2l|, where one could obtain a simultaneous vanishing of the splitting of the axial 
partners and of the topological susceptibility when regularization I was used. In fact, the 
inclusion of strange quarks (N/ = 3) leads to non-linear (mixing) effects since the Ua(1) 
anomaly term is trilinear for three flavors and thus would generate additional mechanisms of 
restoration of axial symmetry, even if a constant anomaly strength is taken into account. In 

1(3 




FIG. 4. Normalized topological susceptibility as function of temperature for several values of the 
anomaly parameter a with regularizations I (left panel) and // (right panel) . 

order to understand the present results, let us remember the comparison of the expressions 
for x within the SU(3) and SU(2) models as discussed in Sec. Ill CI 

A meaningful difference between both models is that in SU(3) we can define an effective 
anomaly strength, <?D e// = 9d (qq), that, as the temperature varies, acquires a temperature 
dependence through the variation with temperature of the quark condensate. Therefore it is 
possible to obtain restoration of chiral axial symmetry without additional assumptions, the 
infinite cutoff at finite T leads to the vanishing of the quark condensate at high temperatures 
and, as a consequence of the effective anomaly strength and of the quantities that depends 
on it: the topological susceptibility and the mass splitting between the axial partners. In 
this case restoration of axial symmetry is a mere consequence of the full restoration of 
chiral symmetry. It is not necessary to implement an anzats related to explicit instanton 
suppression. As we will show, the vanishing of the mass splitting between axial partners in 
SU(2) demands such a mechanism. 

The results plotted in Figs. I2H31 already give an indication that a weaker axial symmetry 
breaking (in this case obtained by hand, by giving low values to a) favors the restoration 
of axial symmetry and reduces the mass splitting of the axial partners (see Fig. [3]). In 
fact, the restoration of this symmetry can only be modeled in a phenomenological way by 
making g 2 temperature dependent. Therefore we will calculate the meson mass spectrum 
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FIG. 5. Left panel: masses of the meson chiral partners (tt, a) and (77, ao) with varying temperature 
for Q2{T) as a decreasing function of temperature, using regularization II (A = const.) at finite T. 
The dotted line represents the qq threshold 2M U . Right panel: normalized mass splitting of axial 
partners as function of temperature for different values of a. 

with gi fixed and #2 a decreasing exponential of the temperature, #2 CO = (?2(0)e~ < - T ~ 300 ^ 20 . 
The result is plotted in Fig. [5] where we can see that, although the effective restoration of 
chiral symmetry (vanishing of the mass splitting between a and tt, and between 77 and ao) 
occurs first, at T = 245 MeV, the degeneracy of the tt, 77, ao and a occurs at T = 285 MeV. 
Therefore, although the analysis of the effects of high momentum quark states and of the 
instanton suppression have similarities in what concerns the characteristic temperatures for 
the phase transitions and the high temperature behavior of the topological susceptibility, 
the evidence that they are indeed mechanisms with a different physical meaning appears 
when we look for a possible degeneracy of the axial chiral partners. 



IV. CONCLUSIONS 

The present study is dedicated to the analysis of the behavior of various observables that 
signal the restoration of the chiral and axial symmetries in regards to temperature, within 
a SU(2) Polyakov-Nambu-Jona-Lasinio model with a 't Hooft interaction term. 

Two types of regularization of finite temperature were considered; the first regularization 
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consists in using the cutoff only on divergent integrals (i) and the other consists in always 
using a finite cutoff (II). The effect of the anomalous coefficient dependent on temperature, 
g2(T), was also analyzed. The restoration of chiral and axial symmetries with temperature 
was studied, considering the effects in the symmetry restoration process of the type of 
regularization applied. A comparative study between the results here obtained and the ones 



calculated within the PNJL model in SU(3) j^l was also performed. Similarly to SU(3), it 
was verified in the present work that the critical temperatures obtained with regularization 
I are closer to lattice results than with regularization II; however, the effective restoration 
of the chiral and axial symmetries does not occur simultaneously. 

We hence conclude that when the coupling constants are considered constant, the restora- 
tion of the axial symmetry, in the SU(2) sector, does not occur as a natural consequence of 
the complete restoration of the chiral symmetry, being this fact a fundamental difference be- 
tween the SU(3) and the SU(2) sectors, where it is not sufficient to apply the regularization I 
for all the observables, related with the U^(l) symmetry, to vanish. It is therefore necessary 
to use an additional mechanism for this restoration, such as a dependence on the temper- 
ature of the coupling coefficient that is related to the anomaly. Once g 2 is the coupling 
constant of the term that breaks the axial symmetry and simulates the instanton effect, this 
result seems to indicate that, in SU(2), the instantons are not completely suppressed by the 
restoration of the chiral symmetry. 

A relevant contribution that the present study offers for the understanding of physics 
associated with the breaking and restoration of the Ua(1) symmetry refers to the fact that 
the analysis of the topological susceptibility, x> * s the "necessary condition", but not the 
"sufficient condition" , for the study and comprehension of the restoration of the axial sym- 
metry. 
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Appendix A 



Some meaningful differences between the results within the SU(3) and the SU(2) models 
appear, so it is useful, in order to understand these differences, to present the basic ingre- 
dients of the SU(3) model. The SU(3) NJL model with the 't Hooft determinant has the 
following Lagrangian: 



C 



a=0 



9D 



det [q(l + j 5 )q] + det [q(l - 75 )g] 
which can be rewritten as: 

C = q^ d ^-m)q+^{(q\ a q)S ab (qX b q) + (q il5 X a q)P ab {qi 1 ,\ b q)) 
where the following projectors have been introduced: 

Sab = 9S$ab + 9DDabc (?A C g) , 

Pab = gs$ab ~ 9DD ahc (q\ c q) ■ 



(Al) 

(A2) 

(A3) 
(A4) 



{qX c q) being the vacuum expectation values. The constants D abc coincide with the SU(3) 
structure constants d a b c for a, b, c = (1, 2, . . . , 8) and D obc = —^5 bc and D 000 — w|. 
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